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CN ' ABSTRACT 

We present the second XMM-Newton observation (85 ks) of the narrow-line Seyfert 1 galaxy 
(NLS1) I Zw 1 and describe its mean spectral and timing characteristics. On average, I Zw 1 
is ~ 35 per cent dimmer in 2005 than in the shorter (20 ks) 2002 observation. Between the 
two epochs the intrinsic absorption column density diminished, but there were also subtle 
changes in the continuum shape. Considering the blurred ionised reflection model, the long- 
On' term changes can be associated with a varying contribution of the power law component 
q relative to the total spectrum. Examination of normalised light curves indicates that the high- 
S_h ■ energy variations are quite structured and that there are delays, but only in some parts of 
^ \ the light curve. Interestingly, a hard X-ray lag first appears during the most-distinct structure 
C3 ■ in the mean light curve, a flux dip ~ 25 ks into the observation. The previously discovered 

broad, ionised Fe Ka line shows significant variations over the course of the 2005 observa- 
tion. The amplitude of the variations is 25 — 45 per cent and they are unlikely due to changes 
in the Fe Ka-producing region, but perhaps arise from orbital motion around the black hole 
or obscuration in the broad iron line-emitting region. The 2002 data are re-examined for vari- 
ability of the Fe Ka line at that epoch. There is evidence of energy and flux variations that are 
associated with a hard X-ray flare that occurred during that observation. 
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1 INTRODUCTION 

I Zwicky 1 (I Zw 1; z = 0.0611) is considered the prototypical 
narrow-line Seyfert 1 galaxy (NLS 1) mainly due to its distinct opti- 
cal properties. At X-ray energies the active galactic nucleus (AGN) 
exhibits behaviour that is somewhat more modest than for typi- 
cal NLSls (e.g. Leighly 1999a, 1999b). However, a short (20 ks) 
XMM-Newton observation conducted in 2002 revealed a plethora of 
interesting characteristics, such as: hard X-ray flaring, unresolved 
transition array (UTA) absorption, and possibly multiple Fe Ka 
emission lines (Gallo et al. 2004a, hereafter G04). 

The discoveries in G04 prompted further examination with a 
second longer XMM-Newton observation. This is the first of three 
papers which describe the X-ray properties of I Zw 1 utilising all 
XMM-Newton data. In this work, we present an 85 ks XMM-Newton 
observation and describe the mean spectral and timing characteris- 
tics of the AGN, as well as compare the data to the earlier 2002 



observation. The complex ionised absorber and its long-term vari- 
ability is discussed in detail in the work of Costantini et al. (2007). 
In addition, I Zw 1 exhibits remarkable short-term spectral vari- 
ability during the 2005 observation. This is presented in Gallo et 
al. (2007). 



2 OBSERVATIONS AND DATA REDUCTION 

I Zw 1 was observed for the second time by XMM-Newton (lansen 
et al. 2001) on 2005 luly 18-19 during revolution 1027 (obsid 
0300470101) (hereafter we will refer to this as the 2005 observa- 
tion). The total duration was 85 ks, during which time all instru- 
ments were functioning normally. The EPIC pn (Striider et al. 2001) 
and MOS (MOS 1 and MOS2; Turner et al. 2001) cameras were op- 
erated in small-window mode with the medium filter in place. The 
Reflection Grating Spectrometers (RGS1 and RGS2; den Herder 
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et al. 2001) also collected data during this time, as did the Optical 
Monitor (OM; Mason et al. 2001). 

The Observation Data Files (ODFs) from the 2002 and 2005 
observations were processed to produce calibrated event lists us- 
ing the XMM-Newton Science Analysis System (SAS v7 . . 0). 
Unwanted hot, dead, or flickering pixels were removed as were 
events due to electronic noise. Event energies were corrected for 
charge- transfer inefficiencies. EPIC response matrices were gen- 
erated using the SAS tasks ARFGEN and RMFGEN. Light curves 
were extracted from these event lists to search for periods of high 
background flaring. Some background flaring is detected toward 
the end of the 2005 observation, and the data during those peri- 
ods have been neglected. The total amount of good exposure time 
selected was 58 ks and 82 ks for the pn and MOS detectors, respec- 
tively. Source photons were extracted from a circular region 35" 
across and centred on the source. The background photons were ex- 
tracted from a larger off-source region and appropriately scaled to 
the source-selection region. Single and double events were selected 
for the pn detector, and single-quadruple events were selected for 
the MOS. The data quality flag was set to zero (i.e. events next to a 
CCD edge or bad pixel were omitted). The total number of source 
counts collected by the pn instrument in the 0.3 — 10 keV range was 
approximately 2.75 x 10 5 . In the 0.5 - 10 keV band, the MOS in- 
struments collected about 8.4 x 10 4 (MOS1) and 8.9 x 10 4 (MOS2) 
source counts. Comparing the source and background spectra we 
found that the spectra are source dominated between 0.2 — 12 keV. 
The EPIC pn spectrum is analysed in the 0.3 — 10 keV band where 
the calibration is most certain. To avoid cross-calibration uncertain- 
ties at low energies, the MOS spectra are neglected below 0.5 keV. 
The 2005 spectra from all instruments show relatively good agree- 
ment in these energy bands, within known calibration uncertainties 
(Kirsch 2006), so they are fitted simultaneously unless stated oth- 
erwise. As the broad-band light curve is less sensitive to calibration 
uncertainties, the data between 0.2 — 12 keV are utilised. 

The RGS were operated in standard Spectro+Q mode in 2005, 
and a total exposure of 85 ks was collected. The RGS data are used 
in the analysis of the ionised absorber in I Zw 1 (Costantini et al. 
2007) and will not be discussed here. 

During 2005, 1 Zw 1 was observed with the OM in Fast-mode 
through four filters (U, UVW1, UVW2 and UVM2). The pur- 
pose was to examine possible optical variability. Unfortunately, the 
AGN was not well centred in the small, fast-mode frame and fell 
close to the detector edge. Consequently, optical variability could 
not be examined in a robust manner. 




Figure 1. The 2005 0.2 - 12kcV EPIC pn light curve in 200 s bins. The 
start time corresponds to the beginning of the pn exposure. 
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3 2005 X-RAY LIGHT CURVES 

During the 20 ks 2002 observation I Zw 1 exhibited a modest X- 
ray flare, which lasted a total of ~ 5 ks (G04). Notably the flare 
appeared to be quite hard, concentrated in the 3 — 12keV band, 
thus inducing spectral variability. Aside from this lone event the 
light curve was relatively quiet in comparison to many NLS 1 . 

The 2005 light curve appears more active, consistently vari- 
able with, on average, larger amplitudes than observed in 2002 
(Figure [TJ. Perhaps the most interesting event in the light curve 
is the sharp drop about 25 ks into the observation. There is spectral 
variability related to the light curve, and this is discussed at length 
in Gallo et al. (2007). Here, we will focus on the flux variations. 

Light curves were created in three EPIC pn energy bands: 
0.2 - 0.7 keV, 2 - 3keV and 3 - 12keV. Using the light 
curves in 1000 s bins, each of the high-energy light curves was 



Figure 2. The CCF from cross-correlating the energy bands shown in each 
panel. In each panel the negative delay curve is reflected about the zero 
axis and shown as a red, dashed curve on the positive side. This is done to 
highlight the asymmetry of the CCFs toward positive lags meaning that the 
soft bands are leading the variations. 

cross-correlated (following Edelson & Krolik 1988) with the the 
0.2—0.7 keV band to search for possible lags (positive delay means 
hard-band lags the soft-band). All of the cross-correlation functions 
(CCFs) were asymmetric toward positive lags (Figure|2j. As will be 
demonstrated in the spectral analysis, the softest band is the region 
where both the power law component and the soft excess compo- 
nent influence the spectrum. Both higher bands are dominated by 
the power law component alone. 

There are a few AGN where time delays as a function of 
Fourier frequency (see Nowak et al. 1999 and references therein 
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Figure 3. Light curves in the 0.2 - 0.7keV (black dots) and 3 - 12kcV 
bands normalised to their respective means. A hard band lag is apparent in 
between ~ 25 — 52 ks (marked by the vertical dashed lines). The start time 
corresponds to the beginning of the pn exposure. 



for details) have been determined (see Papadakis et al. 2001 for the 
first reported detection in an AGN). If there were a single delay be- 
tween the light curves in the various energy bands, then the time 
delay between all of the Fourier frequencies would be the same. In 
contrast, what is typically observed is that as the frequency of the 
Fourier components decreases, the delays increase. Furthermore, as 
the energy band separation between two light curves increases, the 
delays, at the same frequency, also increase (e.g. see Arevalo et al. 
2006 for a recent discussion). 

In the case of I Zw 1, we cannot estimate time lags as a func- 
tion of frequency due to limited signal in the hard-band light curve. 
However, the CCFs are consistent with what we observe in other 
AGN and NLS1. The CCFs show broad humps around zero time 
lag with an asymmetry toward positive lags, but not a single, dis- 
tinct delay. This is consistent with the hypothesis that in I Zw 1, 
just like in other Seyferts, the sine and cosine functions in the hard- 
band light curves are delayed with respect to those in the soft-band 
light curves, but the delay is not the same for all of them: the de- 
lays increase with decreasing frequency and increased separation 
between the energy bands. 

Of interest may be the particularly small values measured in 
the CCFs (0.5 — 0.7). We examine this further by directly com- 
paring light curves in each energy band normalised to their respec- 
tive mean (Figure[3j. Most obvious in the comparison between the 
0.2 — 0.7 keV and 3 — 12 keV band is the significance of the delay 
between ~ 25 — 52 ks (vertical dashed lines in Fig.[3}. In particular, 
the delay seemingly begins at the onset of the sharp dip in the light 
curve. If we calculate the CCF in this time regime alone we see a 
clear dominant lag at ~ 1 ks (Figure|4]l. 

However, as Figure [3] demonstrates, the lag is not always 
present. Sometimes there are no delays, and at other times the light 
curves show different behaviour and even anti-correlations (e.g. 
~ 20 — 25 ks and > 75 ks into the observation). All of these factors 
contribute to the small value of the CCF peaks. Similar behaviour 
was seen in the NLS1 IRAS 13224-3809 (Gallo et al. 2004b) where 
the authors suggested possible alternating lags between high- and 
low-energy bands (i.e. sometimes hard band lags and other times it 
leads). 
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Figure 4. The CCF when cross-correlating the 0.2 — 0.7 keV and 3 — 
12 keV light curves between 25 — 55 ks only. A hard band lag appears 
more dominant in this time band and peaks at ~ 1 ks. 



4 2005 SPECTRAL MODELLING 

The source spectra were grouped such that each bin contained at 
least 20 counts (typically more). Spectral fitting was performed 
using XSPEC vll . 3 . (Arnaud 1996). All parameters are re- 
ported in the rest frame of the source unless specified otherwise. 
The quoted errors on the model parameters correspond to a 90% 
confidence level for one interesting parameter (i.e. a A\ 2 = 2.7 
criterion). A value for the Galactic column density toward I Zw 1 
of 5.07 x 10 20 cm" 2 (Elvis et al. 1989) is adopted in all of the 
spectral fits. K-corrected luminosities are calculated using a Hubble 
constant of Hq=70 km s" 1 Mpc -1 and a standard flat cosmology 
with Om = 0.3 and tt A = 0.7. 

In the analysis of the 2002 data, G04 were able to fit the spec- 
tra with various phenomenological models, but could not distin- 
guish a superior one. While they were able to establish that a sec- 
ond, soft component was a necessary part of the continuum they 
were unable to determine whether the soft excess was curved (e.g. 
blackbody-like) or flat (power law-like). 

4.1 The 2-10 keVspectra 

As was done for the 2002 data (G04), we built-up a broad-band 
spectral model for the 2005 observation utilising all the EPIC data. 
Fitting the 2 — 10 keV EPIC spectra with a single power law (each 
instrument having free photon index (r) and normalisation to ac- 
commodate existing cross-calibration uncertainties) resulted in a 
satisfactory model (x 2 /dof = 1.03/1090). Adding a Gaussian 
profile to this continuum, with a common energy (E) and width 
(a) among all three instruments was a significant improvement 
(A^ 2 = 59.5 for 5 free additional parameters). The best fit line 
parameters were E = 6.70 ± 0.14keV, a = 523+ 2 ™ eV and 
EW pn = 265±*° 5 , EWmi = 220±^°, EWmi = 239±^ eV. 
G04 introduced the possibility of multiple line features in the 
Fe Kq region of the 2002 spectrum, including an intrinsically nar- 
row 6.4 keV line. No such line is present in the 2005 spectrum. A 
narrow, neutral line is completely insignificant in the pn spectrum 
and any such feature in the MOS has an EW < 27 eV. 

A comparable improvement (Ax 2 = 38.0 for 5 free addi- 
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Figure 5. The residuals (data/model) from fitting a model composed of a 
blackbody, power law and Gaussian profile (modified by absorption) to the 
EPIC data (see third row in Table [TJ. Black dots, red open squares, and 
green diamonds correspond to the pn, MOS1, and MOS2 data, respectively. 
The data are rebinned for display purposes. 



tional parameters) can be achieved if a disc line profile (Fabian 
et al. 1989) rather than a Gaussian profile is used for the broad 
residuals. In this case the best-fit energy and inner disc radius are 
E = 6.55^0 ft) keV an d Rin =6 — 26 gravitational radii, respec- 
tively. The disc inclination (i) is fixed at 50°, to remain consistent 
with our finding in Section RTTl Leaving the parameter free to vary 
does not enhance the fit statistics, nor does it affect the best-fit line 
energy. The line strength is EW ~ 200 eV in each spectrum. For 
simplicity, we adopt the Gaussian profile. 



4.2 Simple broad-band models 

Extrapolating the 2-10 keV fit to 0.3 keV in the pn and 0.5 keV 
in the two MOS left a gradual rise in the residuals that is typically 
associated with the soft excess. 

The entire 2005 0.3 — 10 keV band was modelled with various 
continuum models. Through repeated effort we found it necessary 
to include in all models an additional neutral absorber at the redshift 
of the source, as well as a low-energy (~ 650 eV) absorption edge 
that mimics the more complicated ionised absorption. 

As with the 2002 observation, the 2005 data show that a sec- 
ond continuum component is required to fit the spectrum of I Zw 1 
(see Table[TJ. According to an f-test, the addition of a second con- 
tinuum component (i.e. 2 additional parameters) to a single power 
law is highly significant (F-test probability ~ 10 -17 ). Also like 
the 2002 observation, a single superior model does not stand out. 
However, there are similarities among all the models in Table [T] 
One may notice that the second component is not highly curved. 
Even in the case of the blackbody plus power law model, the high- 
temperature of the blackbody component will give it a flatter ap- 
pearance in the EPIC bandpass. 

Of interest is that there is detectable diminishing of the intrin- 
sic column density in I Zw 1 since 2002. The exact change depends 
on what continuum model is assumed, but appears to be at least on 
the 50 per cent level. 

The residuals from the blackbody plus power law fit are shown 
in Figure[5] Notable in Figure[5]are the positive residuals, which re- 



main in the fit above ~ 8 keV. These residuals are apparent in all 
broad-band fits to the 2005 data (see also Figure [7] and the lower 
panel of Figure [8]l , with the exception of the double power law 
model. The most likely cause is the significant spectral variability 
exhibited by I Zw 1 in 2005 (Gallo et al. 2007). Gallo et al. suggest 
that the power law photon index could be changing during the 2005 
observation, which would explain the difficulty in fitting the aver- 
age spectrum with a single power law. The reason that the residuals 
are minimal in the double power law model could simply be that 
the second power law does not begin to dominate the spectrum un- 
til rather high energies (~ 5.5 keV in the pn spectrum). 

Based on the blackbody plus power law model, the average 
0.3 — lOkeV flux during 2005, corrected for Galactic absorption, 
is ~ 1.4 x 10 -11 erg cm -2 s _1 . The luminosities in the 0.3 — 
10 keV and 2-10 keV bands are 1.28 x 10 44 erg s _1 and 0.43 x 
10 44 erg s _1 , respectively. I Zw 1 is ~ 35 per cent dimmer than it 
was in 2002. For comparison, the estimated 2 — 10 keV luminosity 
during the ASCA observation was 0.53 x 10 44 ergs' 1 (Leighly 
1999b). 



5 A COMPARISON OF THE 2002 AND 2005 SPECTRA 
5.1 Summary of the 2002 observation 

The 0.2 — lOkeV spectrum of I Zw 1 has long been known to 
differ from those of other NLS1. It possesses a rather weak soft 
excess and is modified by intrinsic neutral and ionised absorption. 
The 2002 observation of I Zw 1 (G04) was rather typical of what 
had been seen previously (e.g. Leighly 1999a, b). 

There were several simple continuum models (e.g. black- 
body plus power law, broken power laws, multiple blackbodies 
plus power law) that when modified by intrinsic ionised absorp- 
tion (an absorption-edge), neutral absorption (absorption column 
density) and an Fe Kq emission line, fit the EPIC spectra ade- 
quately in 2002. None of these continuum models stood out as su- 
perior. Since the 2002 data were reprocessed with the newest cal- 
ibration we refit the simple models to the newly processed data 
and find comparable results. The one modification was in the in- 
trinsic absorption column density, which was slightly lower with 
the new calibration (A^h = 4.5 ± 0.5 x 10 20 cm -2 as opposed to 
Nh = 10.0 ± 0.1 x 10 20 cm -2 ). 

A particularly interesting spectral result from the 2002 obser- 
vation was the possible presence of multiple iron lines. In addition 
to a broad and ionised line, the spectrum also required a narrow, un- 
resolved emission line at 6.4 keV . The most obvious interpretation 
adopted by G04 was that this feature was produced in distant mate- 
rial such as the torus, but an alternative interpretation is offered in 
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5.2 Spectral changes since 2002 

The 2005 spectrum does appear substantially different from the 
2002 observation. The broad-band differences between 2002 and 
2005 can be depicted by extrapolating the best-fit 2 — 10 keV power 
law to both spectra down to lower energies (Figure[6j. 

The 2005 observation shows a much softer spectrum, a sug- 
gestion that the level of absorption may have changed since 2002. 
The best-fitting power law in the 2 — 10 keV band is slightly flat- 
ter in 2005 (r = 2.15 ± 0.02 compared to V = 2.28 ± 0.03 in 
2002), which likely also contributes to the apparently stronger soft 
excess at that epoch. The inset in Figure[6]shows the residuals in the 
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Table 1. Fit parameters for the simple models (column 1) applied to the 2005 full-band data. Fit quality is given in column 2. Neutral absorption intrinsic to 
I Zw 1 is given in column 3 ( X 10 20 cm -2 ). In column 4, the continuum parameters are shown. If the parameter is not linked between the various instruments 
it is identified as pn, Ml (MOS1) and M2 (MOS2). E b is the break energy (in keV) for the broken power law model. The photon index above (r^) and 
below (r s ) the break energy are also shown. The blackbody temperature (kT) is given in eV. Column 5 provides the absorption-edge parameters: energy (E) 
and optical depth (r). All fits except for the partial covering model (second row) include an FeKa emission line with comparable parameters to the results 
presented in the text. The line was not statistically necessary in the partial covering model. The parameters of the partial covering absorber, column density 
(JVupc, x 10 22 cm -2 ) and covering fraction (/), are given in column 6. 
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Figure 6. Top panel: Best-fit power law modelled to the 2002 (red marks, 
T = 2.28±0.03) and 2005 (black marks, T = 2.15±0.02) 2-10keVpn 
spectrum, and extrapolated to 0.3 keV. Lower panel: Resulting ratio from 
the above model. Top panel (inset): Ratio in the 4 — 9.5 keV band, binned 
for display purposes. 

Fe Kq band. As with previous observations, the iron line appears 
broad and ionised. The agreement in the appearance of the iron line 
at both epochs is good. 



We applied the blackbody plus power law model from 2002 to 
the 2005 data to examine consistency between the two epochs. Sim- 
ply rescaling the model was unacceptable (x„/dof = 5.7/846). 
Allowing the level of neutral absorption to vary freely significantly 
improves the fit quality (xt/doi = 1.27/845) indicating that the 
column density of the neutral absorber has dropped in 2005 to 
Nn ~ 1.4 x 10 20 cm -2 . Allowing the low-energy absorption edge 
to vary in energy and depth is a further improvement (A^ 2 = 36.2 
for 2 additional free parameters). The overall fit is not acceptable 
(X^/dof = 1.23/843) and demonstrates that subtle changes in the 
intrinsic continuum shape since 2002 may be appropriate. 

If rather we consider a constant absorber and allow only the 
continuum parameters to vary, we obtain a worse fit (xl/doi = 
1.35/843). Indeed variable absorbers seem necessary and are sup- 
ported by the RGS data (Costantini et al. 2007). 



6 A REFLECTION INTERPRETATION FOR I Zw 1 
6.1 The 2005 observation 

The spectrum of I Zw 1 is power law dominated, but there is the 
necessity for a small contribution from a soft component below 
~ 2 keV. The ionised disc reflection model (Ross & Fabian 2005) 
blurred by relativistic effects has had relatively good success in 
duplicating the appearance of the soft excess (e.g. Crummy et al. 
2006). 

All the EPIC spectra were fitted with the blurred reflection 
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Figure 7. Top panel: The best fit blurred reflection model (only the pn 
model is shown for clarity). The reflection (red, dash-dot curve) and power 
law (blue, dashed curve) components, which make up the combined fit 
(black, solid curve) are presented. Bottom panel: The residuals of the model 
to the EPIC spectra. Symbols are as defined in Figure[5] 



Table 2. The blurred reflection and power law model (with absorption) ap- 
plied to the 2005 EPIC spectra (x^/dof = 1.10/1627). The absorption 
parameters are as defined in Table[T] The blurring parameters are the emis- 
sivity index (a); inner {Rin) and outer (Rout) disc radius in gravitational 
radii (r g = GM/c 2 ); and disc inclination (i). The iron abundance and 
ionisation parameter (f; = AizFn -1 , where F is the total incident flux and 
n is the hydrogen number density) of the reflector are also given. 
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plus power law model. As in the previous fits, we included an in- 
trinsic neutral absorption component and an absorption edge as a 
pseudo-ionised absorber. Since the model is complex, to simplify 
the fitting process the normalisations of the MOS spectra were 
linked to the pn spectrum and only a constant scale factor was used 
to correct for possible flux offsets between instruments. The inci- 
dent spectrum for the reflection component was linked to that of the 
power law component. 

The fit quality is not as satisfactory as some of the simple mod- 
els (x^/dof = 1.10/1627), but it does have the advantage of of- 
fering a physical explanation for the origin of the soft excess and 
broad Fe Ka line. Fit parameters are shown in Table [2] and the 



Figure 8. Residuals from a self-consistent power law plus blurred reflection 
model fitted simultaneously to the pn data from 2002 (top panel) and 2005 
(bottom panel). 



model along with fit residuals are presented in Figure|7] The frac- 
tion of reflected-to-total flux in the 0.3 — lOkeV range (corrected 
for Galactic and source absorption) is ~ 0.30. Drawing conclu- 
sions from the fit parameters should probably be done with cau- 
tion given the complexity of the low-energy absorption. However, 
it can be noted that the disc inner radius is not extremely small 
{Rin — 5.0+^'g r g ) and the emissivity profile is not very steep 
(a — 2.481J5 gi)- Both points suggest that the reflection component 
does not necessarily originate from the innermost confines around 
the black hole as is suggested for many NLS1. 

The advantage of this model is that it can account for the iron 
emission and soft excess in a self-consistent manner. Moreover, as 
will be discussed in Gallo et al. (2007), this reflection scenario can 
potentially explain some of the complicated spectral variability ex- 
hibited by I Zw 1 . 



6.2 Simultaneous fit to the 2002 and 2005 data 



Based on the investigation in Section [5721 most of the spectral dif- 
ferences between 2002 and 2005 can be associated with changes 
in the level of absorption. However, there do appear to be subtle 
changes in the shape of the continuum as well. As I Zw 1 was 
nearly twice as bright in 2002, in terms of the blurred reflection 
model the changes in the continuum shape could be interpreted as 
differing contributions from the direct power law component. 

The data from both epochs were fitted simultaneously to ex- 
amine consistency in the reflection scenario. For simplicity, many 
parameters were fixed to the best-fit 2005 values (see Table [3j, 
The primary parameters that were left free for each dataset were 
the power law and absorption parameters, and the normalisation of 
the reflection component. The iron abundance and disc inclination 
were linked between the two spectra as large changes in these pa- 
rameters would not be expected over 3-years. A reasonable fit can 
be achieved (xt/doi — 1.06/1480) in this manner (see Figure [8] 
and Tabled. 

As is seen in Table[3]fhere is no strong indication for a change 
in the shape of the direct power law component between 2002 and 
2005. In addition, most of the blurring and reflection parameters 
have been fixed. The most notable changes are in the amount of 
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Table 3. The blurred reflection and power law model (with absorption) 
applied simultaneously to the 2002 and 2005 pn spectra (xi I d°f = 
1.06/1480). The parameters are as defined in Tableland / marks the ones 
that are fixed in the fit. The 2002 and 2005 values are given in columns (2) 
and (3), respectively. The reflection (R) and flux fractions (Fr) (as defined 
in the text) are given in the final two rows. 



(1) 

Model Component 


(2) 
2002 


(3) 
2005 


Absorption 
(neutral and 
ionised) 


N m = 4.38l°;2° 

E = 623 ± 10 
, _ o o +0.03 


0.22^0;^ x 10 20 cm- 2 
GMtl 1 eV 
0.23 ± 0.02 


Power law 


T = 2.28 ±0.02 


n n, +0.10 

Z ' Z± -0.08 


Blurring 


a = 2A8 f 
R in = 5.0-f 
Rout = 200' 
i = 58 ± 3 


2.48* 
5.0* r g 

58 ±3 


Reflector 


Fe= 1.10 ±0.13 

£ = soeo-f 


1.10 ±0.13 solar 
3060-^ erg cms -1 


R 
Fr 


0.27 ±0.04 
0.22 ±0.03 


0.45 ± 0.04 
0.31 ± 0.03 



absorption (as noted in previous fits and by Costantini et al. 2007), 
the reflection fraction (R, fraction of reflected emission to power 
law emission in the 0.3—10 keV band), and in the level at which the 
reflection component contributes to the total 0.3—10 keV spectrum 
(F R = 0.22 ± 0.03 in 2002 and 0.31 ± 0.03 in 2005). 



7 IRON LINE VARIABILITY 

7.1 FeKa variability during the 2005 observation 

In the analysis of G04 there were indications that the broad feature 
in the 2002 observation was composed of multiple (at least two) 
narrower emission lines. In the 2005 observation, the residuals can 
be fit by a single Gaussian profile. Based on the width of the line 
(a — 523^57 eV) in the EPIC spectra, the line-emitting region is 
between 12 — 66 r s . Similar distance estimates come from the disc 
line fit and the blurred reflection model. Given the long observa- 
tion and proximity of the line emitting region makes it plausible to 
search for line variability. 

The mean pn spectrum was divided into 17 segments, each 
of 5 ks. The exposure in each segment was ~ 3.5 ks, taking ac- 
count of the live-time of the pn camera in small window mode. The 
average number of source counts in the 5.0 — 7.2 keV band was 
slightly more than 300. The counts were grouped into spectra such 
that each bin contained at least 20 counts. For each 5 ks spectrum 
the continuum was determined by fitting the 4—10 keV band with 
a power law, while excluding the 5 — 7.2 keV range. The omitted 
band was replaced and the spectrum was rebinned into channels 
of 250 eV. Excess counts above the continuum in each 250 eV bin 
were recorded. From this information an excess residual image was 
created as illustrated in Iwasawa et al. (2004; and in prep). The re- 
sult is analogous to modelling each spectrum with a power law plus 
Gaussian profile and measuring the flux in the Gaussian compo- 
nent. The advantage of creating an excess residual image is that we 
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Figure 9. Excess-residual image for the 2005 pn spectrum in the 4.5 — 
8.2 keV band. Bright spots indicate positive residuals in the spectrum, 
above a power law continuum. The time resolution is 5 ks, and the spec- 
tra are binned into 250 eV wide channels. The bright band between 5.8 — 
6.7 keV corresponds to the broad emission line seen in the mean spectrum. 
The feature appears to be variable over the course of the observation. 




20000 40000 60000 80000 

Time (s) 

Figure 10. The light curve of the broad Fe Ka emission in I Zw 1 during 
the 2005 observation. Bins are 5ks. The l-<r uncertainties were estimated 
as described in the text. A constant fit is shown for comparison. 

can circumvent the limitations of working in the Poisson regime, 
thus making trends in the variations of the line-flux more visible. 

The resulting image from the 2005 observation is shown in 
Figure [9] and clearly indicates that the broad-line flux is variable. 
Most notably, the line flux diminishes significantly 45 — 60 ks from 
the start of the observation. 

To create a light curve from the image, the excess counts be- 
tween 5.8 — 6.7 keV were accumulated. Uncertainties on the light- 
curve data points were estimated from Monte Carlo simulations of 
the excess residual image. One-thousand images were created in an 
identical manner as done for I Zw 1. The spectra from which the 
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Figure 11. Top panel: The 6 — 10 koV light curve of the continuum used to 
establish the iron line variations, normalised to its average flux. The FeKa 
data shown in Figure [Tol are normalised to the average FeKa flux and over- 
plotted (red, dashed curve). Bins are 5 ks. The 1-cr uncertainties on the con- 
tinuum data points are taken from the 6 — 10 keV count rate. Lower panel: 
The CCF of the two light curves shown in the top panel. A positive delay 
indicates that the continuum changes lead the emission line variations. 



images were derived were simulated for identical conditions as the 
2005 observation, assuming a power law plus stationary Gaussian 
profile. The uncertainty in the observed light curve was taken to be 
the standard deviation in the excess counts resulting from the simu- 
lated images. Details of the process including caveats are discussed 
in Iwasawa et al. (2004). 

The Fe Ka light curve shows clear variability (FigurefTUt with 
a constant fit being unacceptable (xt — 2.45). The fractional vari- 
ability amplitude (F V3bI , using uncertainties derived in Edelson et 
al. 2002) is 34.8 ± 9.7 per cent. 

A light curve of the 6 — lOkeV continuum was generated 
based on the power law fits used to derive the Fe Kq light curve. 
The uncertainties in the 6—10 keV flux were assumed to be the 
same as the uncertainty in the corresponding count rate. The con- 
tinuum light curve, normalised to its average flux, is shown in Fig- 
ureQT]with the broad-line flux light curve (Figure [Tot, normalised 
to its average flux, overplotted as a dashed curve. Interestingly, the 



10000 20000 

Time (s) 




(Ph/3600s/cm**2) 

Figure 12. Excess-residual image for the 2002 pn spectrum in the 4.5 — 
8.2 keV band. Bright spots indicate positive residuals in the spectrum, 
above a power law continuum. The time resolution is 3.6 ks, and the 
spectra are binned into 250 eV wide channels. The bright band between 
6.0 — 6.5 keV corresponds to the broad emission line seen in the mean 
spectrum. Between 10 — 15 ks there is also brightening in the 5.5 — 6.0keV 
band. 

curves appear to be anti-correlated in some parts of the light curve. 
Cross-correlating the broad-line and 6 — 10 keV continuum light 
curves indicates a weak correlation with the broad-line light curve 
lagging the continuum fluctuations by ~ 15 ks (lower panel Fig- 
ure [TTJ. 

The fractional variability in the 6 — lOkeV continuum is 
Fvar = 20.8 ± 4.0 per cent, less than the F var in the Fe Ka line. 
Statistically, the amplitude of the continuum variations is consis- 
tent with the Fe Kq variations within about 2-cr. However, the dif- 
ference is likely larger since the uncertainties on F var are conser- 
vatively over-estimated (Edelson et al. 2002). 

7.2 Revisiting the Fe Ka variability in 2002 

In fitting the spectra in various segments of the 2002 observa- 
tion, G04 found little evidence of intensity variations in the broad- 
line. We take advantage of the excess-residual technique to re- 
examine the broad-line variability in the 2002 observation. An 
excess-residual image has been created for the 2002 pn observa- 
tion in a nearly identical manner as was done for the 2005 obser- 
vation. The only difference is that the time resolution is slightly 
better during 2002 because the mean spectrum was divided into 6 
segments of 3.6 ks each (rather than 5ks). This can be done with- 
out degrading the data quality because the I Zw 1 was brighter in 
2002 and because the observation was conducted in large-window 
mode, which has a CCD live-time of ~ 95 per cent; therefore the 
actual exposure in each segment of the 2002 observation is almost 
the same as the exposure in each segment of the 2005 data. 

The 2002 excess-residual image is shown in Figure [T2] The 
broad-line is dominant in the 6.0 — 6.5 keV range, but there is also 
apparent brightening that occurs in the 5.5 — 6.0 keV band between 
10 — 15 ks. This is of particular interest because the continuum light 
curve shows a hard X-ray flare that peaks during this time interval 
(see G04). 
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Figure 13. Various light curves from the 2002 observation of I Zw 1 
normalised to their respective averages. The black, dashed curve is the 
3 — 12keV continuum light curve highlighting the hard X-ray flare re- 
ported previously by G04 (uncertainties in the curve are found in G04). The 
filled, red circles connected by a solid line mark the light curve extracted 
from the residual image (Figure [T2l in the 5.5 — 6.0 keV range (i.e. the red 
line). The peak intensity in the red line corresponds with the peak intensity 
in the 3 — 12 kcV continuum light curve. The open, blue squares connected 
by a dotted line is the light curve extracted from the residual image in the 
6.0 — 6.5 keV band (i.e. the blue line). 

Light curves were extracted from the excess-residual image in 
the 6.0 — 6.5 keV (blue line) and 5.5 — 6.0 keV (red line) bands. 
They were then normalised to their respective means and plot- 
ted with the normalised 2002 3 — 12 keV continuum light curve 
(Figure Q/J}. As the observation was so short uncertainties in the 
residual-image light curves were not derived in a robust manner as 
was done for the 2005 data. We also note that propagation of Pois- 
son uncertainties in each data bin yields very large error bars (on 
the order of 100 per cent), therefore we are only illustrating a trend 
in Figure [T3l 

The light curves suggest that the variations are more intense in 
the 5.5 — 6.0 keV-line than in the 6.0 — 6.5 keV-line. There also is 
no obvious connection between the 6.0 — 6.5 keV-line light curve 
and that of the continuum. On the other hand, the 5.5 — 6.0 keV 
variations appear to be well correlated with the high-energy contin- 
uum variations and, despite the unknown size of the true error bars 
in the 5.5 — 6.0 keV-line light curve, the trend is obvious. 



8 DISCUSSION 

8.1 I Zw 1: a not so prototypical NLS1 in X-rays? 

The strong Fe II emission in the optical spectrum of I Zw 1 was 
noted a half century ago (Sargent 1968), and given the intensity 
and narrowness of the lines the AGN became a laboratory for Fe II 
studies (e.g. Phillips 1976, 1977, 1978; Boroson & Green 1992; 
Laor et al. 1997; Vestergaard & Wilkes 2001). Its narrow Balmer 
lines were also noted fairly early on (e.g. Oke & Lauer 1979). 

When a formal definition of the NLS1 class was established 
(Osterbrock & Pogge 1985), based on optical spectral properties, 
it was realised that I Zw 1 was already the best-studied member of 
this class of AGN, likely leading to its labelling as the prototype 
NLS1. However, some X-ray properties are also characteristic of 



NLSls and I Zw 1 does not exhibit all these qualities to noticeable 
levels. In fact, in some instances I Zw 1 even displays behaviour 
that is atypical of the class. For example, I Zw 1 is observed through 
a column density intrinsic to the AGN. This is unusual for NLSls 
(e.g. Boiler, Brandt & Fink 1996), but also for type 1 AGN in gen- 
eral. 

The soft excess emission above the underlying steep power 
law is weak in I Zw 1 compared to most NLSls, and not because of 
the relatively high intrinsic absorption (e.g. Leighly 1999b, Tanaka 
et al. 2005). Even the clear detection of a broad Fe Kq emission 
line in I Zw 1 seems unusual given that NLSls appear not to show 
broad emission lines so obviously (e.g. Fabian et al. 2002, 2004; 
Guainazzi et al. 2006). 

The identification of I Zw 1 as the prototype NLS 1 originated 
due to its well-studied optical properties, but it now lingers mainly 
as a historic perspective. It can even be argued that some of its X- 
ray properties are atypical of the norm. 

8.2 Time lags 

The average CCFs calculated in Section 3 show a tendency for pos- 
itive lags (hard band following soft band). This phenomenon has 
been observed in other AGN (e.g. Arevalo et al. 2006 and refer- 
ences therein) and Galactic black holes (e.g. McClintock & Remil- 
lard 2005). However the origin of the lags is not certain. 

In a simple corona with a uniform temperature and optical 
depth, variations can result from fluctuations in the intensity of the 
seed photons. Variations of this nature should produce fixed delays 
between low- and high-energy bands, resulting in a common time 
delay for all Fourier components. We do not estimate time lags as 
a function of frequency in I Zw 1 ; however the broad, asymmetric 
CCFs suggest that we do not have a single, well defined delay be- 
tween the various energy bands. Examination of normalised light 
curves in different energy bands further supports this. A time lag 
is only seen in parts of the light curves. In fact, at times the light 
curves are not correlated, or are actually anti-correlated. 

Lags are seen in Cyg X-l, which can be attributed to spec- 
tral evolution in magnetic flares as they detach from the accretion 
disc (e.g. Poutanen & Fabian 1999). Another possibility is that the 
lags are manifested by material that propagates toward the black 
hole from regions of soft X-ray emission to regions of hard X-ray 
emission (e.g. Bottcher & Liang 1999; Kotov et al. 2001). 

Based on the observation of another NLS1, 
IRAS 13224+3809, Gallo et al. (2004b) claimed a lag that 
alternated (sometimes hard band followed other times it led), with 
possible dependency on the average flux of the source. There is 
likely similar behaviour here as the first notable lag coincides with 
a sharp flux dip in the light curve about 25 ks into the observation. 



8.3 Diminished neutral absorption since 2002 

Significant variations in the low-energy spectra of some Seyfert 1 
galaxies have been previously observed and interpreted as changes 
in the low-energy continuum component (e.g. Page et al. 2002, 
Starling et al. 2004). In the case of I Zw 1, this does not appear 
to be the most straightforward explanation. There is detectable di- 
minishing of the intrinsic column density in I Zw 1 since 2002. The 
exact change depends on what continuum model is assumed, but 
appears to be at least on the 50 per cent level. 

Changes in the neutral column density have been reported in 
several type 2 (i.e. absorbed) AGN (e.g. Risaliti et al. 2002), but to 
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the best of our knowledge not in a NLS 1 . I Zw 1 is a merging and 
starburst system so there are several sources of absorption that exist 
in the galaxy and could feasibly change (e.g. due to bulk motion 
or dissipation of an absorbing component) over the span of a few 
years, which separates the two observations. 

8.4 Origin of the Fe Kq variability 

8.4.1 Flux variations in the 2005 observation 

The broad Fe Ka line in I Zw 1 shows remarkable short-term vari- 
ability. Understanding the origin of this variability is challenging. 

Figure QT] suggests that the emission line and continuum vari- 
ations could be correlated if we consider the iron-line fluctuations 
to lag the continuum by 10 — 12 ks. Such a delay would place the 
iron-line emitting region at approximately 180 r s in I Zw 1. The 
distance is too large to account for the width of the iron feature 
which based on various spectral fits (e.g. Gaussian profile, disc line 
profile, blurred ionised reflector), originates within about 60 r g . A 
second problem is that the amplitude of the variations appears to 
be larger in the Fe Ka line than in the 6—10 keV continuum that 
would be producing it. Even if we assume that the amplitudes of 
the variations are identical, it would be unlikely that all of the in- 
cident continuum will be reprocessed by the iron producing region 
and then emitted back into our line-of-sight. 

Another possibility is that the variations are a result of some 
extrinsic effect rather than intrinsic changes in the line-emitting re- 
gion. For example, motion within the accretion disc could account 
for some of the flux changes (e.g. Iwasawa et al. 2004). In this 
case, one would expect energy modulation in the excess residual 
map, which are not present in Figure [9] However, the strength of 
the modulations is dependent on the distance of the line-emitting 
region from the black hole and the disc inclination. It could be that 
in I Zw 1 this combination of parameters does not generate observ- 
able modulations. 

In a speculative approach, we fit the Fe Kq line light curve 
(Figure [Tot with a sine curve. The fit was statistically good (xt = 
1.02) and the implied period was 49.1 ks. If the 49.1 ks corre- 
sponds to the orbital period of the emitting region around the black 
hole, then we can estimate the distance (e.g. following Bardeen et 
al. 1972) to be ~ 25 r g . This distance is consistent with that esti- 
mated from the width of the broad Fe Ka line in the spectral mod- 
els, but it is not clear why only the line-emitting region is affected 
in this manner and not the continuum. 

Finally, the Fe Kq variations could be due to variations in 
some optically-thick absorber in the line-emitting region. Such ob- 
scuring material could be created from instabilities in the accretion 
disc of highly inhomogeneous flows (e.g. Merloni et al. 2006). 

8.4.2 Energy and flux variations in the 2002 observation 

The most impressive feature in the 2002 continuum light curve of 
I Zw 1 was a modest X-ray flare concentrated in the 3 — 12 keV 
band. The 6.0 — 6.5 keV line-emission is not particularly variable 
and demonstrates no connection to the continuum variations. This 
behaviour is typical of other results (e.g. Iwasawa et al. 2004). 
However, the 5.5 — 6.0 keV line-emission exhibits more intense 
variations which mimic the continuum light curve. In particular, 
the peak intensity in the 5.5 — 6.0 keV-line light curve corresponds 
well with the peak of the hard X-ray flare. 

Due to limited data, examination of time-resolved spectra dur- 
ing 2002 do not yield significant results. The likely correlation be- 



tween the continuum and 5.5 — 6.0keV-line light curve suggests 
that this is emission due to Fe Ka, which is produced in a spot on 
the accretion disc that is illuminated by the X-ray flare (i.e. perhaps 
co-rotating with the disc). The 6.0 — 6.5 keV-line emission is repre- 
sentative of the average broad-line emission, which is produced at 
distances of ~ 35 r g as suggested by its width. The 5.5 — 6.0 keV- 
line emission and hard X-ray flare are generated closer to the black 
hole where Doppler and gravitational effects are more intense than 
the region where the 6.0 — 6.5 keV emission line is produced (i.e. 
< 35 r g ). 

It is worth noting that the low-energy emission could be re- 
lated to the narrow 6.4 keV (6keV observed-frame) emission fea- 
ture that was identified in the 2002 spectrum by G04. Since the 
feature was not detected in the 2005 spectrum its origin could be 
connected to the transient 5.5 — 6.0keV-line (perhaps the blue- 
wing of this line), which is produced during the X-ray flare. 

The combination of a hard X-ray flare, and redshifted and 
broadened Fe Kq emission line has been seen in a very simi- 
lar object, Nab 0205+024 (Gallo et al. 2004c). However, due to 
limited data, time-resolved spectral modelling did not yield in- 
formation about possible line variability. Re-examination of the 
Nab 0205+024 data making use of excess residual imaging may 
be fruitful. 

8.5 Physical models for I Zw 1 

The level of absorption (neutral and ionised) appears to have 
changed between the two epochs. However, even after consider- 
ing this there does seem to be some subtle change in the intrinsic 
continuum shape. 

It is of interest to compare the changes in the context of a 
physical model such as reflection and light bending. In terms of the 
blurred reflection model, most of the intrinsic changes can be un- 
derstood as arising from different relative contributions of the two 
continuum components. At both epochs (and historically), I Zw 1 
appears to be in a power law dominated state - the soft excess is 
notably weak in this object. 

If the power law component is a compact isotropic emit- 
ting source located above the black hole, as depicted in the light- 
bending scenario (e.g. Miniutti & Fabian 2004), then the changes 
in the continuum shape are largely motivated by the proximity of 
this source to the black hole. In the higher flux 2002 observation, 
the contribution from the reflection component was less, and hence 
the soft excess was weaker. The AGN was brighter and in a power 
law dominated state. In 2005, 1 Zw 1 was dimmer and the reflection 
component (i.e. the soft excess) was stronger. 

This is almost certainly a simplified description of the true sit- 
uation. One would expect other spectral changes to occur as the 
power law component varies its distance from the black hole. For 
example, with increasing distance the photon index should harden 
(i.e. if due to Comptonisation), ionisation parameter may decrease 
as the accretion disc is less illuminated, and Ri n will probably in- 
crease. 

As the 0.3 — 10 keV spectrum of I Zw 1 is dominated by one 
component and subject to ionised absorption, it would not be sur- 
prising if other physical models provided adequate fits. Broad-band 
absorption models, like the blurred absorption models (e.g. Gierlin- 
ski & Done 2004), have not been examined here, but could likely 
duplicate the spectrum of I Zw 1. It may be possible to break this 
degeneracy by examining the spectral variability, which is complex 
in I Zw 1 and is somewhat unlike what is seen in other NLS 1 (Gallo 
et al. 2007). 
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9 CONCLUSIONS 

In this work, we presented an 85 ks XMM-Newton observation and 
described the mean spectral and timing characteristics of the AGN, 
as well as compared the data to a previous, shorter XMM-Newton 
observation from 2002. The complex ionised absorber is discussed 
in detail in the work of Costantini et al. (2007). The remarkable 
spectral variability is presented in Gallo et al. (2007). The main 
findings from this work are the following: 

(1) I Zw 1 is ~ 35 per cent dimmer in 2005 than it was in 2002. 
Between the two epochs the level of neutral column density intrin- 
sic to I Zw 1 changed, as well as the ionised absorber. In addition, 
subtle changes in the continuum shape are apparent. If modelled 
with the blurred reflection model, the long-term continuum changes 
can be understood as arising from different relative contributions of 
the two continuum components (i.e. direct power law and ionised 
reflection). 

(2) Variability is prevalent in all X-ray bands. The most distinct 
structure in the mean light curve is a flux dip ~ 25 ks into the ob- 
servation. Average CCFs show a tendency for positive lags mean- 
ing soft band changes lead harder bands. However, examination of 
normalised light curves indicate that the high-energy variations are 
quite structured and lags are only apparent in some parts of the light 
curve. Interestingly, the hard X-ray lag first appears at the onset of 
the flux dip in the light curve. This possible alternating lag has been 
proposed for another NLS 1 (IRAS 1 3224-3809, Gallo et al. 2004b) 
and it is not consistent with Comptonisation from a single corona 
with a uniform temperature and optical depth, and variable seed 
photons. 

(3) The broad ionised Fe Ka line shows significant variations over 
the course of the 2005 observation, which are not clearly due to 
intrinsic changes in the Fe Ka-producing region. Other, viewer- 
perspective effects, such as orbital motion or absorption in the 
broad iron line-emitting region, could be invoked to describe the 
variability. 

(4) The broad-line variability during 2002 is re-examined here. 
There is evidence suggesting that a relativistically broadened 
Fe Ka line is produced as a result of a hard X-ray flare illuminating 
the inner regions (i.e. < 35 r g ) of the accretion disc. 
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